arXiv: 1503.00625v2 [nucl-th] 21 Apr 2015 


Toward relativistic mean-field description of N-nucleus reactions 


T. Gaitanos M. Kaskulov ^ 

^ Department of Theoretical Physics, Aristotle University of Thessaloniki, GR-54124 Thessaloniki, Greece 
^ Institutfur Theoretische Physik, Universitdt Giessen, D-35392 Giessen, Germany 

email: tgaitano@auth.gr 


Abstract 

In this work we study the antinueleon-nueleus optieal potential in the framework of the non-linear 
derivative (NLD) model with momentum dependent mean-fields. We apply the NLD model to 
interaetion of antinueleons (N) in nuelear matter and, in partieular, to antiproton seattering on 
nuelei. In nuelear matter a strong suppression of the N-optieal potential at rest and at high 
kinetie energies is found and eaused by the momentum dependence of relativistic mean-fields. 
The NLD results are consistent with known empirical N-nucleus observations and agree well 
with antiproton-nucleus scattering data. This makes the NLD approach compatible with both, 
nucleon and antinucleon Dirac phenomenologies. Furthermore, in nuclear matter an effective 
mass splitting between nucleons and antinueleons is predicted. 
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1. Introduction 

An interpretation of recent astrophysical observations on compact neutron stars [1,2] has 
driven the nuclear physics and astrophysics communities to detailed investigations of the nu¬ 
clear equation of state (EoS) under conditions far beyond the ordinary matter [3]. Theoretical 
and experimental studies on heavy-ion collisions [4, 5, 6, 7, 8, 9, 10] over the last few decades 
concluded a softening of the EoS at high densities, which was compatible with phenomenolog¬ 
ical [11], as well as with microscopic models [12]. However, the recent observations of pulsars 
with masses of 1.97 ± 0.04 Mq [1] and 2.01 ± 0.04 Mq [2] gave some controversial insights 
on the high-density EoS. These observations provide the upper limit for the neutron star mass by 
excluding a soft EoS at high baryon densities. 

On the other hand, the in-medium antinucleon-nucleon interactions are closely related to 
the high baryon density domain of the EoS. In highly compressed and cold matter the baryons 
can reach very high Eermi momenta. Such kinematical conditions can be also achieved in an 
antiproton-nucleus reaction for a short time before annihilation. In fact, in antiproton-nucleus 
collisions, where in rare events antiprotons penetrate deep inside the nucleus, significant com- 
pressional effects are expected [13]. Therefore, the study of antinucleon-interactions in nuclear 
medium can provide useful information and additional constraints on the nuclear EoS. 
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The traditional approach for the description of nuclear matter is a well established relativistic 
mean-field (RMF) theory [14, 15, 16, 17, 18]. However, as shown in Refs. [19, 20, 21], RMF 
cannot properly describe in-medium antinucleon interactions using the G-parity arguments only, 
which is supposed to be a symmetry of RMF. On the contrary, the microscopic models for the 
description of in-medium antinucleon potentials strongly rely on G-parity in construction of the 
bare antinucleon-nucleon interaction starting from the nucleon-nucleon potential [22,23]. This is 
considered as a problem in RMF approach, therefore, different solutions have been proposed. For 
instance, in Ref. [20] a violation of G-parity in nuclear matter has been considered. In Ref. [24] 
another method was invented to resolve this issue by using energy dependent interactions. How¬ 
ever, this approach has limitations in density and lacks the thermodynamic consistency. 

In this work we address the issue why the conventional RMF models are not consistent with 
antiproton-nucleus Dirac phenomenology. Our studies are based on the non-linear derivative 
(NLD) model to RMF as formulated in Ref. [25]. The NLD model describes simultaneously 
the density dependence of the nuclear EoS, the momentum dependence of the proton-nucleus 
Schrddinger-equivalent optical potential and the recent observations of neutron star masses. In 
this work we extend our studies from Ref. [25] and show that in a G-parity conserving NLD 
framework the in-medium proton and antiproton optical potentials can be consistently repro¬ 
duced in agreement with empirical data. We further compare the NLD calculations with scat¬ 
tering data. It is demonstrated, in particular, that the NLD model can describe the antiproton- 
nucleus reaction cross section data fairly well. Furthermore, we discuss other effects related to 
the behavior of antinucleons in nuclear matter and predict a considerable in-medium splitting 
between nucleon and antinucleon effective masses. The relevance of our results for upcoming 
experiments at FAIR is discussed. 

2. The NLD Model 

At first, we briefly discuss the NLD model and for additional details we refer to Refs. [25, 26, 
27]. The NLD approach is based on a field theoretical formalism of relativistic hadrodynamics 
(RHD). The NLD Lagrangian 
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contains the free Lagrangians for the Dirac field 4/ and for the a, oj and p meson fields. The 
essential difference to the conventional RHD shows up in the NLD interaction Lagrangian [25] 
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where pm stands for the scalar-isoscalar a, vector-isoscalar and vector-isovector me¬ 
son fields with m = {a, uj, p}, with obvious notations for couplings pm and vertices F^ = 
{ 1, 7 ^, (r are the isospin matrices). 
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In Eq. (3) P, V are non-linear derivative operators, see Ref. [25]. They generate (infinite) 
series of higher-order derivative terms in the generalized Euler-Lagrange equations of motion. 
This series of terms lead to a momentum (energy) eut-off, whieh regulates the high momentum 
behavior of relativistie mean-fields in the interaetion Eagrangian. The important feature of NED 
model is that, independent of the funetional form of the non-linear operators, all the higher-order 
derivative interactions can be resummed exactly in the RME approximation. 

The in-medium RME interactions are obtained from the equations of motion for all degrees 
of freedom. Eor the generalized functional in Eq. (3) the field-theoretical formalism has been 
developed in [25] and is applied in this work. In the RME approximation the Dirac equation for 
nucleons {Ni) forming the nuclear matter reads 

MiS" - E') - (m - E,)| <I>„, = 0. (4) 

Eor antinucleons (iVj) interacting with nuclear matter the Dirac equation is derived from the 
requirement of G-parity invariance 

- E'‘) - (m - E,)] = 0. (5) 

The nucleon and antinucleon vector self-energies are given by 

-f gpTiP^V , (6) 

^vi = + QpTip^V (7) 

where now Ti = +1 for protons (i = p) and r* = — 1 for neutrons (i = n). Note the opposite signs 
in the isoscalar-vector interactions in Eqs. (6) and (7) between nucleons and antinucleons and 
that both, the vector and scalar components of the self-energy, depends explicitly on the particle 
momentum p through the NED regulator function V = V{p) in momentum representation. The 
scalar self-energy in Eqs. (4) and (5) reads as 

= g^aV . ( 8 ) 

Though the scalar self-energies between nucleons and antinucleons are formally the same, the 
kinetic energy dependence will generate a difference between them. Therefore, we distinguish 
them by different symbols in the following. 

The single particle energies E for nucleons and antinucleons are obtained from the corre¬ 
sponding in-medium mass-shell conditions 

En. = 77i*2 ^ (9) 

= + (10) 

with the in-medium (effective or Dirac) masses of nucleons and antinucleons, given by m* = 
m — Tjg and m* = m — respectively. They depend explicitly on particle momentum. 
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In nuclear matter the NLD equations of motion for the meson-fields are redueed to standard 
RMF algebraic equations 
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where the various quantities and Lorentz-densities are given in Ref. [25]. The densities on the 
r.h.s in above equations depend explieitly on the partiele momentum. The meson-field equations 
show a similar form as in standard RMF, however, the essential differenee between the NLD 
approaeh and eonventional RMF appears in the various Lorentz-sourees of the meson fields, 
whieh get regulated at high momenta. 

The way one regulates the high-momentum behavior of fields is phenomenologieal. In 
Ref. [25] a simple monopole form 


A2 - 1 - ^2 ’ 


(14) 


has been used where A is a cut-off parameter. It has been found that Eq. (14) was able to accom¬ 
modate the saturation properties of nuelear matter and, at the same time, the energy dependenee 
of the Schrodinger-equivalent optieal potential. Note, that in the limit A —)■ cx) one has "D —1, 
such that the original RHD Lagrangian is reeovered. Meanwhile in Refs. [32, 33] the NLD model 
has been tested with different choices of the eut-off function. 

The NLD approaeh is thermodynamieally consistent [25]. The parameters of the model (cou¬ 
pling constants and cut-off) were adjusted to the empirieal properties of nuelear matter and to the 
momentum dependenee of the in-medium potential at saturation. As discussed in detail in [25], 
the NLD approaeh ean be reliably extrapolated to very high densities. As an example, in Lig. 1 we 
demonstrate our new ealeulations of the equation of state (EoS) (pressure versus energy density, 
main panel) for nuclear matter in /)-equilibrium. The NLD ealeulation (solid eurve) reproduces 
the empirical analyses from astrophysieal studies (shaded areas) fairly well [3, 28]. Of partieular 
interest are the constraints from neutron star measurements. This is shown in the inner part of 
Lig. 1 in terms of the mass-radius relation of neutron stars, where the NLD result is eompared 
with various astrophysieal observations (double NS systems and J1903-1-0327), and in particular, 
with the recent measurements on binary pulsars [1,2]. 

The explieit momentum-dependence of the NLD mean-fields is essential for a simultaneous 
description of low and high density observations involving, in particular, the empirical results 
from Dirae-phenomenology on the Sehrodinger-equivalent optieal potential. The latter is impor¬ 
tant for the extension of our studies to in-medium antinueleon interaetions. 
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Figure 1: (Main panel) Pressure P as function of the energy density e in the NLD model. The shaded areas are 
the results of empirical studies [3, 28]. (Inserted panel) Mass-radius relation for neutron stars as calculated in NLD 
model, together with results from astrophysical observations and constraints; the four horizontal shaded bands refer 
to measurements from double neutron star (NS) systems [29, 30], from the pulsars PSR J1903H-0327 [31], PSR 
J1614-2230 [1] and PSR J0348H-0432 [2]. The other shaded areas bordered by thick curves indicate parameter space 
excluded by general relativity, causality (shaded area on the top-left) and rotational constraints (shaded area on the 
bottom-right) [29, 30]. 


3. Results and discussion 

The parameters of the NLD model used here have already been fixed in Ref. [25] at the 
bulk properties of nuelear matter and in-medium proton optieal potential. Therefore, we do 
not readjust any parameters in this work when extending the NLD model to the deseription of 
antinueleons in nuelear matter. We just use the G-parity arguments for the eonstruetion of the 
in-medium antinueleon interaetion. 

We start the diseussion with the in-medium antiproton sealar and veetor self-energies. In 
Fig. 2 they are shown as a funetion of the in-medium antiproton kinetie energy 

Ekin = Ep-m= - m (15) 

at three baryon densities (note that is negative for antinueleons by negleeting the isovee- 
tor eontributions, see Eq. (7)). The eonventional Waleeka approaeh (NL3 parametrization [34], 
dashed eurves) shows the usual linear behavior in density for the veetor self-energy, while the 
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Figure 2: Kinetic energy dependence of the scalar and vector Lorentz-components of the antiproton self-energy in 
nuclear matter at densities of p = psat (left), p = 2psat (middle) and p = 3psat (right) using the NL3 model (dashed 
lines), the NL3 model with rescaled couplings with the factor 0.25 [20] (dash-dotted) and the NLD approach (solid 
lines). 


scalar component saturates with density. Obviously, the Waleeka self-energies are flat as a func¬ 
tion of the kinetic energy. On the other hand, the NLD relativistic self-energies (solid curves) 
deerease with inereasing kinetie energy at eaeh baryon density. The strong suppression of the 
high-density and high-momentum tails of the NLD self-energies result from the momentum de- 
pendenee of the fields, see Eqs. (6,7). In partieular, at zero kinetie energy the Waleeka model 
leads to a value of — S 5 —700 MeV,whioh remains eonstant in energy and is very deep. On 
the contrary, the NLD model reduces considerably the deepness of the in-medium antinueleon 
potential in all ranges of shown kinetie energies. 

As found in Ref. [20], in order to reproduee the data from antiproton-indueed reactions within 
a transport model, the antinueleon-meson eouplings of the Waleeka model had to be resealed by a 
phenomenologieal faetor of 0.2 — 0.3. In Fig. 2 we show sueh ealeulations in the Waleeka model 
with resealed eouplings (faetor 0.25, dash-dotted eurves in Fig. 2). As one ean see, the resealed 
Waleeka model [20], whieh is supposed to explain the antiproton-nueleus data, just reproduces 
in average the NLD results. 

The momentum dependenee of the relativistie mean-fields ean be probed by looking at the 
Sehrodinger-equivalent optieal potential. Its real part is given by 

fi^ef/opt = -V-S+^{S^-V^) , (16) 

m 2m 

where E = (Ej^.) is the energy of a nueleon (antinueleon) with bare mass m inside nuclear 

matter at a fixed baryon density, V = ST (^J is the eorresponding veetor self-energy and 
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Figure 3: Real part of the Schrodinger-equivalent optical potentials for antiprotons (main panel) and protons (insert) 
interacting with nuclear matter at saturation density as function of the kinetic energy. The shaded areas in the main 
panel and the filled circles in the inserted are extractions from empirical analyses [20, 35, 36]. The dashed (dot- 
dashed) curves are calculations in the (rescaled) Walecka model and the solid curves are calculations within the 
NLD approach. 


S = Eg (Ss) denotes the sealar self-energy eomponent. 

Fig. 3 shows 93et/opt as funetion of the kinetie energy for in-medium antiprotons at saturation 
density. The shaded areas for antiprotons around zero kinetie energy and antiprotons at high 
energies are the results of phenomenologieal analyses of Refs. [36, 20]. It ean be seen that the 
momentum dependenee is erueial for a eorreet deseription of the deepness of the optieal potential 
(solid eurve), whereas the standard Waleeka approaeh (dashed eurve) does not reproduee the 
absolute value and the trend of empirieal data. The eomparison with the empirieal analyses is 
improved in the resealed Waleeka model (dot-dashed eurve). 

The inner graph in Fig. 3 shows the eorresponding results for the proton ease, where the NLD 
ealeulations (solid eurve) deseribe the Dirae phenomenology very well. However, the Waleeka 
model (dashed eurve) again disagrees with the data (filled symbols [35]). We stress again that 
the NLD model parameters are kept fixed when going from in-medium proton to in-medium 
antiproton interaetions. 

At this point we would like to explain the behavior of the optieal potential in the NLD model 
at kinetie energies around i?kin — 0 MeV. This ean be done by eonsidering the in-medium 
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dispersion relation, that is 


m = Vp" + (17) 

where p = |p | and is equal to 

p" = -(S,)"+(E^p)" + 2m(E,-E;p) (18) 

with Ss (^p) being the sealar (veetor) RMF self-energies of the antiproton at fixed nuclear matter 
density. Hence, neglecting isovector contributions the single-antiparticle momentum in nuclear 
matter contains the sum of the scalar and the vector components of the RMF self-energies, p ~ 
Ss — S„p. Note, that in the case of in-medium protons the sign of the vector self-energy changes, 
p ~ Ss — S°p. In the antiproton case the single-antiparticle momentum is very big, while in 
the proton case it is essentially given by the small difference between scalar and vector mean- 
fields. This is demonstrated in Fig. 4, where the in-medium kinetic energy E — m is shown 
as a function of the single-(anti)particle momentum for nuclear matter at saturation density of 
Psat = 0.155 fm“^. Indeed, while in the proton case the momentum for a particle with in¬ 
medium kinetic energy Ekin = 0 MeV is of the order of the Fermi momentum only, its value 
increases largely for the antiproton case. For an antiproton in matter at saturation density it 
approaches values of several hundreds of MeV, namely p{Ehin = 0) ~ 750 MeV, which are 
comparable with the nucleon mass. 

This momentum-shift shown in Fig. 4 has significant consequences for the behavior of the 
potential deepness at zero in-medium kinetic energy, iRtUoptiEkin = 0), when the self-energies 
depend on single-particle momentum. In fact, as it is shown in the inserted graph in Fig. 4, 
the NLD self-energies are reduced with increasing single-particle momenta. At a momentum of 
p ~ 750 MeV the values of = 163 MeV and S^p = 134 MeV for the scalar S and vector 
V self-energies, respectively, are obtained. The sum of + S^p = 297 MeV is considerably 
smaller than the corresponding sum at Fermi momentum. Since with increasing momentum the 
NLD fields are getting reduced, and in the antinucleon-case the sum of scalar and vector self¬ 
energies matters, therefore these effects results in the reduced value of the in-medium antinucleon 
potential. 

Another effect generated by the momentum shift between in-medium protons and in-medium 
antiprotons, concerns the effective particle masses. Fig. 5 shows the effective masses for protons, 
m* = m — Ss, and antiprotons, rn* = m — in nuclear matter versus their kinetic energies 
and at two baryon densities. A mass shift between protons and antiprotons is observed, which 
becomes pronounced with decreasing kinetic energy and increasing baryon density. At the same 
value of the kinetic energy the proton scalar self-energy is obtained at a smaller proton mo¬ 
mentum than in the antiproton case, where is larger. Therefore, at the same kinetic energy the 
effective proton mass is smaller as compared to the effective antiproton mass. This mass split¬ 
ting between protons and antiprotons becomes moderate with increasing kinetic energy because 
of the less pronounced momentum dependence of S* at high particle momenta. However, the 
mass splitting is strongly density dependent. This is demonstrated in Fig. 6, where the effec¬ 
tive proton and antiproton masses are shown at zero kinetic energy as a function of the baryon 
density. 


8 



1000 



0 250 500 750 1000 1250 

p [MeV/c] 

Figure 4: In-medium kinetic energy of a proton (dashed curve, NLD p) and an antiproton (solid curve, NLD p) as 
function of single-particle momentum for nuclear matter at saturation density psat = 0.155 fm~^. The calculations 
are performed in the NLD model. The inserted panel shows the momentum dependence of the scalar (solid curve) 
and vector (dashed curve) self-energies of a proton in nuclear matter at psat- 


One would expect that the only effect of the G-parity is the sign change in the vector self¬ 
energy. This is indeed the case in the standard Walecka models where there is no difference 
between the in-medium proton and antiproton masses. However, the NLD self-energies depend 
on momentum, as supported by Dirac phenomenologies. The impact of momentum dependence 
supplemented by different momentum shifts of protons and antiprotons in the scalar self-energies 
produces this observed mass splitting. Note the analogy with the isospin induced mass splitting 
between protons and neutrons in isospin-asymmetric nuclear matter. In the NLD model there is 
a mass splitting effect between nucleons and antinucleons also in G-asymmetric nuclear matter. 

Next, we compare the NLD calculations with scattering data. This requires the knowledge 
of both, the real and the imaginary part of the in-medium optical potential. Using the dispersion 
relation one can calculate the imaginary part of Uopt from the real part in Eq. (16). The imaginary 
part of the Schrddinger equivalent optical potential is given by [37] 


JmUopt(p) = -—V 

71 


^tUoptip' 


rpf 2 _ 


-dp', 


(19) 


where p = \p \ stands for the antiparticle momentum and V denotes the principal value. The 
real part (Hef/opt is taken from the NLD model. Note that Eq. (19) makes sense only if (Hef/opt 
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Figure 5: Effective mass of a proton (dashed curves) and an antiproton (solid curves) as function of the kinetic 
energy in nuclear matter at ps = Psat (left) and at pB = ipsat (right) in the NED model. 


does not diverge with inereasing momentum like in the standard Waleeka model. This property 
of 9fe(7opt is realized in the NLD model. 

The imaginary part of the in-medium antiproton optieal potential is shown in Fig. 7 as func¬ 
tion of the antiproton kinetic energy. The empirical data are taken from the analyses on stopped 
antiprotons, see for details Refs. [38]. The NLD calculations are in agreement with the empirical 
analyses within their uncertainties. The depth of the imaginary part of the optical potential is 
also consistent with phenomenological studies of antiprotonic atoms [38]. 

Having both, the real and imaginary parts of the in-medium potential as function of beam 
momentum and density, we can now use the NLD optical potential to estimate the antiproton 
total cross sections in the scattering off nuclei. The numerical procedure is as follows. For a 
given nuclear target we use the empirical Wood-Saxon density distributions [43]. They provide 
us with the density profiles which are used as an input to calculate the corresponding nuclear 
potential (7opt as a function of the radial coordinate. Then, at intermediate energies the scattering 
amplitude / is obtained using the the eikonal approximation, that is 

poo 

f{^) = ip dbbJo{Ab)il - . (20) 

Jo 

The momentum transfer is given by A = 2p sin id/2, Jq denotes the Bessel function, r = b± + zn 
and the eikonal x reads 

X = Uopt{b,z)dz. (21) 
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Figure 6: Effective mass for a proton (dashed curve) and for an antiproton (solid curve) at i?kin = 0 MeV in nuclear 
matter as a function of the baryon density in the NED model. 


Using the optical theorem the total cross section takes the form 

47r 

atot = —= 0 ), ( 22 ) 

P 

where the imaginary part of the amplitude is given by 

poo 

Jm/(t))=p / (i66Jo(A6) (l — cos . (23) 

Jo 

In Fig. 8 we show the NLD calculations for the total cross section as function of the antiproton 
beam momentum on different nuclear targets. The total antiproton cross sections are obtained 
from experimental absorption cross section data including optical model analyses for the missing 
contributions to the total one. In particular, the data symbols for the nuclear targets C, Cu and 
Al at low beam momenta are taken from Ref. [39], whereas the experimental cross sections at 
the higher energies for C and Cu are from [40]. For antiproton reactions on the heavier Ph 
target shown in Fig. 8 there exist only few experimental estimations [41, 42]. As one can see, 
the calculations reproduce the absolute values, momentum and mass number dependence of the 
total cross sections fairly well. 

The comparison of the NLD model with the results of empirical studies and with exper¬ 
imental data is remarkable. Our results compare well also with calculations done within the 
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Figure 7: Imaginary part of the Schrodinger-equivalent optical potential (see Eq. (19)) as function of antiproton 
kinetic energy at saturation density. NED calculations (solid curve) is compared with empirical extrapolations 
(symbol, vertical lines) taken from [38]. 


Giessen-BUU transport model [20]. In that work the imaginary part of the in-medium antipro¬ 
ton interaetion was determined with the help of the eollision integral of the transport equation, 
where the elastie, inelastie and annihilation processes were taken into account. There a standard 
RMF model has been also used. However, as we have already discussed, in Ref. [20] one had 
to strongly reduce the antinucleon-meson coupling constants to achieve a good description of 
antiproton absorption cross section data at low and intermediate energies. Therefore, in line with 
Ref. [20] a possible violation of G-parity was discussed. This conclusion is not supported by the 
results of the present work, where a simple constraint imposed by the energy dependence of the 
Schrodinger-equivalent optical potential may result in a consistent G-parity conserved descrip¬ 
tion of both, the in-medium nucleon and antinucleon interactions. 

4. Summary 

The NLD model with momentum dependent mean-fields has been applied in the past suc¬ 
cessfully for the descriptions of nuclear matter at low and high baryon densities. In this work the 
NLD formalism has been extended to the in-medium interactions of antinucleons. We studied 
the optical potential of antiprotons inside nuclear matter and then used it for the description of 
antiproton reactions off nuclei. 
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Figure 8: Total antiproton reaction cross sections as function of beam momentum for various target nuclei, as 
indicated. The NLD calculations (solid curves) are compared with experimental data (symbols), taken from [39, 40, 
41, 42] (see text for more details). 


It has been demonstrated that momentum-dependent interactions are able to describe simul¬ 
taneously the nucleon and antinucleon optical potentials in agreement with available empirical 
data. In particular, we find a considerable reduction of the optical potential for antiprotons as 
compared to the standard RMF models. The same momentum dependence results in a new ef¬ 
fective mass splitting between nucleons and antinucleons in nuclear matter. This sort of medium 
induced G-asymmetry becomes pronounced for low energy antinucleons and increases with ris¬ 
ing baryon density. 

Furthermore, the imaginary part of the in-medium optical potential has been determined by 
using the dispersion relation. The results concerning the imaginary part of the antiproton poten¬ 
tial are compatible with available phenomenological analyses. 

We then applied the NLD optical potential to the description of antiproton-induced reactions 
on nuclear targets. Our calculations reproduce fairly well the dependence of the total reaction 
cross sections on beam momentum and nuclear mass number. 

We demonstrated that momentum dependent interactions, as realized in NLD model, lead 
to a consistent description of empirical information and experimental data. Together with the 
successful application of NLD approach to low and high density nuclear matter, we conclude 
that, a proper relativistic mean-field formalism should account for the momentum dependence of 
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fields. This is important for a correct RMF description of nuclear matter and finite systems. 

As an outlook, it would be interesting to apply the NLD in reactions with heavy-ion and an¬ 
tiproton beams within transport theoretical descriptions. In particular, the mass splitting between 
nucleons and antinucleons, as predicted in this work, is expected to show up in various observ¬ 
ables. For instance, it is known that in- and out-of-plane collective flows as well as subthreshold 
particle production are sensitive to possible mass splitting effects. Thus, the investigations of 
compressed matter in reactions relevant for PANDA at FAIR may provide interesting signals as 
probes of the NLD predictions. 
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